t 


i 


Technical  Report  ixixjxxxXx 


A  RANDOM  WALK  TREATMENT  OF 
NEUTRON  DIFFUSION  IN  SLABS 


September  1966 


NAVAL  FACILITIES  ENGINEERING  COMMAND 


U.  S.  NAVAL  CIVIL  ENGINEERING  LABORATORY 
Port  Hueneme,  California 


£00¥0?d-£>0j.i/ 

Dufribofion  of  'his  docum.nt  i,  unlimit.d 

D  D  Yr=r^ 
V^nnnv?m 

OCT  18 1966  m 

,  u  J 


t 


A  RANDOM  WALK  TREATMENT  OF  NEUTRON  DIFFUSION  IN  SLABS 

1L .  , 

Tiecljnical  Report  R-482  ^  ^ 


Y-R01 1-01-01-101 


L.  Eaton  and 


C. 


M. 


Huddleston 


i  ** 


mT  & 

tt 


.'it 


ABSTRACT 

The  concept  of  a  Markov  chain  has  been  used  to  treat  the  random  processes  of 
scatter  and  absorption  which  occur  when  neutrons  are  incident  on  a  slab  of  finite 
thickness,  such  as  a  shield  or  an  inside  wall  of  a  shelter  entranceway.  It  is  assumed 
that  scattering  is  isotropic  in  the  laboratory  system  and  that  the  scattering  and  absorp¬ 
tion  cross  sections  do  not  change  during  a  neutron-scattering  history.  The  encouraging 
results  obtained  to  date  indicate  that  for  thin  slabs  the  analytical  random  walk  method 
may  have  important  advantages  over  Monte  Carlo  calculations  (which  require  lengthy 
computer  runs  to  obtain  acceptably  small  statistical  variances),  moments  method 
calculations  (which  are  actually  appropriate  only  for  infinite  medium  cases),  and 
numerical  solutions  of  the  neutron  transport  equation  (which  are  lengthy,  tedious, 
and  necessarily  approximate). 


The  Laboratory  invites  comment  on  this  raport,  particularly  on  tha 
results  obtained  by  those  who  have  applied  the  information. 
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INTRODUCTION 

In  calculating  shielding  effectiveness  it  is  often  desired  to  have  values  for 
neutron  albedo.  In  particular,  the  following  items  are  of  interest  for  a  neutron 
incident  on  a  slab  of  isotropically  scattering  material: 

1.  Its  total  albedo;  i.e.,  its  probability  of  eventually  backscattering  out  of 
the  slab. 

2.  An  analysis  of  the  total  albedo  into  differential  albedo  as  a  function  of 
the  polar  angle  of  return  or  backscatter  direction. 

3.  An  analysis  of  the  total  into  the  integer  numbers  of  impacts  experienced 
before  returning. 

4.  The  probability  density  of  the  depth  in  the  slab  at  which  a  backscattered 
neutron  makes  its  last  impact. 

The  classical  partial  differential  transport  equation  uses  time  as  an  independent 
variable.  It  appeared  that  a  possible  approach  to  the  foregoing  set  of  problems 
would  be  to  write  and  solve  a  psuedo  transport  equation  in  which  the  time  scale  is 
replaced  by  a  scale  linear  in  the  number  of  historical  neutron  impacts.  Difficulties 
in  writing  such  an  equation  were  discussed  with  Dr.  Klaus  H.  Daniel,  Professor, 
Department  of  Mathematics,  University  of  Maryland,  during  his  stay  at  the  U.  S. 
Naval  Civil  Engineering  Laboratory  in  the  summer  of  1965.  Dr.  Daniel  suggested 
the  possibility  of  obtaining  insight  by  first  obtaining  some  approximate  solutions 
using  increments  of  slab  thickness  (instead  of  differential  thickness)  and  then  studying 
these  approximate  solutions  in  the  hope  of  learning  how  to  write  the  correct  equation 
to  go  to  the  limit  in  the  distance  variable  normal  to  the  slab  face. 

This  suggestion  was  pursued.  There  is  still  no  guarantee  that  the  psuedo 
equation  will  be  written;  nevertheless,  the  method  for  obtaining  approximate  solutions 
by  dividing  slab  thickness  into  increments  seems  to  hold  promise  as  a  usable  method 
in  its  own  right.  The  description  of  this  method  constitutes  the  purpose  of  this  report. 


STATEMENT  OF  THE  PROBLEM 

Consider  the  case  of  a  neutron  incident  at  some  polar  angle  to  the  surface  of  a 
plane  slab  of  material.  The  slab  has  a  finite  thickness  and  consists  of  matter  which  can 
absorb  neutrons  or  scatter  neutrons  isotropically  in  the  laboratory  system.  The  probabil¬ 
ities  for  absorption  and  scatter  at  each  collision  are  constant;  i.e.,  the  energy  of  the 
neutron  does  not  change.  The  foregoing  description  is  equivalent  to  the  case  of  a 
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broad,  parallel  beam  of  thermal  neutrons  incident  on  a  plane  slab  of  a  moderator  whi 
has  some  absorption.  It  is  required  to  find  a  probabilistic  description  of  the  behavior 
of  the  neutrons  in  the  slab. 

Mathematically,  the  problem  can  be  phrased  in  terms  of  those  quantities  which 
are  given  and  those  quantities  which  are  desired. 

Given: 

1.  A  plane  slab,  of  thickness  T,  has  known  values  for  E^,  1^,  and  I  ,  which 

are,  respectively,  the  total  cross  section,  the  scattering  cross  section,  and  the  absorp' 

tion  cross  section.  L  =  E  +  I 
t  s  a 

2.  Scattering  is  isotropic  in  the  laboratory  system. 

3.  The  slab's  dimensions  are  such  that  it  occupies  the  space  (-ao  <  X  <  ao), 

(-00  <  Y  <  ao),  (0  ^  Z  *  T).  The  face  at  Z  =  T  will  be  named  the  forward  face,  and 
the  initial,  near,  or  backward  face  is  at  Z  =  0. 

4.  One  neutron  is  directed  from  the  region  (Z  <  0)  toward  the  origin  (on  the 
backward  slab  face  at  Z  =  0)  in  a  course  which  makes  a  given  angle  0Q  with  the 
-Z  axis,  and  an  angle  (tr/2  -  0Q)  with  the  -X  axis. 


1.  Total  albedo,  (?  (probability  of  eventually  scattering  back  to  Z  <  0). 

2.  Differential  albedo,  a,  as  a  function  of  0,  so  that 

2n  ir/2 

I  I  a(0)sin0  d0d#  =  A 

0  0 

3.  Analysis  of  total  albedo  into  the  number  of  impacts  experienced;  i.e.. 


w 

’  I  4 


4.  Analysis  of  total  albedo  into  depth,  Z,  (0  s  Z  £  T),  probability  density  of 
the  depth  at  the  point  of  last  impact  before  returning. 

5.  Total  absorption,  Q  (probability  of  even'ual  absorption  of  the  neutron  by 
the  slab). 

6.  Analysis  of  fl  into  total  number  of  impacts  experienced;  i.e., 


CD 

■  l  a° 
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7.  Analysis  of  Cf°  into  depth  probability  density  of  point  of  absorption. 

8.  Probability,  (?*  ,  that  the  neutron  will  go  through  the  slab  eventually. 

9.  Analysis  of  ff'  into  the  number  of  impacts  experienced;  i.e. , 


a' 


where  (?„  is  the  probability  of  uncollided  transmission. 

10.  Analysis  of  ( (t  -  (1^)  into  depth  density  of  last  impact  before  leaving 
forward  face  at  Z  =  T . 

11.  Expected  number  of  neutron  impacts,  E(N). 

12.  Expected  number  of  impacts,  E(N^),  if  neutron  is  backscattered  from 
the  slab. 

13.  Expected  number  of  impacts,  E(N0),  if  neutron  is  absorbed. 

14.  Expected  number  of  impacts,  E(N$),  if  neutron  is  scattered  through  the  slab. 

15.  Expected  number  of  impacts,  E(Nj.),  if  eventually  neutron  goes  through 

the  slab. 


For  the  sake  of  clarity,  a  few  obvious  consequences  of  the  mathematical  model 
are  now  stated. 

Granted  that  an  impact  is  definitely  in  the  offing,  the  probability  that  the 
neutron  will  be  absorbed  at  that  impact  is  2^,/Ij .  Granted  that  it  does  not  become 
absorbed  at  that  impact,  its  scattering  direction  density  is  uniformly  ]/4n  in  all 
directions,  so  that 


sphere 


or,  equivalently. 


I 


sin©  d0d<£ 
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a  +  r  +  a'  =  i 


3 


E(N)  =  Q  E(Nb)  +  (f  E(N0)  +  ( Cf '  -  Qx0  )E(NS) 

=  a  E(Nb)  +  a°  E(N0)  +  a*  E(Nt) 

oo 

-  E  jp(j> 

j-o 

OG 

E(N>  =  y;  jp<j) 
j='i 

A  PRELIMINARY  RESULT 

Given  a  sphere  with  infinite  radius  and  of  material  with  total  cross  section 
1^,  given  that  isotropic  scattering  of  a  nonabsorbed  neutron  has  just  occurred  at  a 
known  impact  point,  and  given  that  an  arbitrarily  directed  W  axis  with  origin  at 
that  known  impact  point  has  been  erected;  then,  both  the  scattering  angles  and  the 
distance  to  the  next  impact  point  are  random  variables.  Also,  the  W  coordinate 
(perpendicular  projection  onto  the  W  axis)  of  the  next  impact  point  is  a  random 
variable  over  the  range  (-ao  <  W  <  ao).  Let  the  units  for  reckoning  distance  be  the 
same  in  the  W  system  as  in  the  (X,  Y,  Z)  system. 

Assuming  the  foregoing,  consider  the  following  problem:  Find  the  density 
function,  h(W),  for  the  W  coordinate  of  the  next  impact.  For  motivation,  it  is 
mentioned  now  that  in  sections  to  follow  the  W  axis  will  be  oriented  parallel  to 
the  Z  axis  for  the  slab  problem. 

Because  this  whole  report  depends  upon  the  validity  of  this  function,  h(W) 
will  be  sought  by  two  different  methods:  first,  directly;  and  second,  indirectly, 
by  computing  first  the  distribution  function  H(W),  and  then  differentiating  it  with 
respect  to  W  to  obtain  h(W). 

Before  starting  either  solution,  some  preliminary  comments  seem  appropriate. 
By  a  distribution  function,  H(W),  is  meant  a  function  bounded  below  by  zero  and 
above  by  unity  such  that 

W 

H(W)  =  j*  h($)d£ 

-co 

and  dH(W)/dW  =  h(W),  everywhere  that  H(W)  is  differentiable,  where  h(W)  is  the 
density  function  of  the  random  variable  W. 
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Because  of  symmetry  about  the  point  W  =  0,  H(0)  =  1/2,  and  h(W)  =  h(-W). 
Finally,  the  following  general  properties  of  all  distribution  functions  are  noted 

H(-ao)  =  0 

H(ao)  =  1 

H(W)  is  monotonically  nondecreasing  with  increasing  W. 

Direct  Derivation  of  h(W) 


Consider  any  ray  from  the  origin  W  =  0,  the  immediately  previous  impact 
point,  into  the  right  half  infinite  sphere  W  >  0,  making  an  angle  0  <  n/ 2  with  the 
positive  W  axis.  For  a  neutron  traveling  this  ray,  the  W  coordinate  of  the  next 
impact  is  a  random  variable  with  density  function  g(W)  =  Zj,  sec  0  e”^  Wsec0^ 

W  >  0.  Since  this  is  true  for  aJJ  such  rays  in  the  right  half  sphere,  their  uniformly 
(1/4 ff)  weighted  "sum"  is 


2p  Y2  -Ij  Wsec  0 

h(W)  =  J  J  sec  0  e  sin  0  d0  drj 

0  0 


where  7]  represents  angles  in  the  plane  perpendicular  to  the  W  axis  through  the 
point  W  =  0.  This  recognizes  the  fact  that  scatters  from  this  impact  point  into  the 
left  half  sphere  or  the  right  half  sphere  are  equally  likely. 

Thence,  by  integrating  over  7J, 

Z  Y2  -Z  Wsec0 

h(W)  =  y  J  tan0e  d0 

0 


To  facilitate  this  remaining  integration,  the  following  change  of  variable  is  made: 
V  =  Z^  W  sec  0 

dV  =  Z^Wtan0sec0  d0 
dV  _  , 

-y-  “  tan  0  d0 
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Substitution  gives 


oo 

Z.  -V 


hW  =  y  |  * 


h(W)  =  -  y  Ei(-ZfW)  W  >  0 


By  symmetry. 


h(W)  =  -  y  Ei{ZfW)  W  <  0 
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Values  for  the  exponential  integral  Ei(>)  are  tabulated.1' 


Indirect  Derivation  of  h(W) 
Noting  that 


H(0)  = 


and  that  for  any  neutron  leaving  the  origin  W  =  0  and  entering  the  right  half  sphere, 
the  probability  of  its  impacting  again,  before  its  projection  on  the  W  axis  becomes 
as  large  as  any  particular  positive  W  value,  is  1  -  e”H  Wsecp^ 

i  \"c  ,  -I,W!ec(\ 

H(W)  =  -j  +  f  J  (l  -  .  )sin0d0 

0 

*{?  -Z  WseciS 

=  1  -  y  J  e  sin£d/3  W  >  0 

0 


=  1  - 


NL  p  -V 

^  I  y 

ztw  V 
-*tw  zw 


=  1  - 


y-  Ei(-ZfW)  W  >  0 
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INTRODUCTION  TO  MARKOV  CHAIN  CONCEPTS 

Consider  some  phenomenon  in  which  a  "particle"  at  any  one  instant  may  be  in 
any  one  of  a  known  number  (few,  several,  many,  or  a  countably  infinite  number)  of 
different  possible  states  (positions  or  conditions).  Given  that  the  particle  can  change 
its  state  sequentially  in  a  chain  of  steps,  and  that  for  each  existing  state,  the  proba¬ 
bility  of  its  being  in  any  other  state  immediately  after  the  next  step  is  known  for  all 
states;  then  the  problem  of  computing  for  a  particle  in  one  known  given  state  the 
probability  of  its  being  in  some  particular  given  state  exactly  n  steps  later  is  a  Markov 
chain  problem. 

The  following  example  may  clarify  the  foregoing  abstraction.  Bill  has  four 
coins  and  Joe  has  three  coins,  and  they  agree  to  gamble  these,  one  at  a  time  (say  by 
matching  or  flipping)  where  for  each  person,  the  probability  of  winning  or  losing  the 
coin  each  time  (each  step)  is  one-half.  From  Joe's  point  of  view  there  are  eight 
possible  states:  broke,  has  one  coin,  ...,  has  six  coins,  has  seven  coins.  The  chain 
of  steps  is  the  successive  encounters  or  matchings  with  the  attendant  change  in  wealth. 
It  is  mathematical!/  convenient  to  indicate  the  different  possible  states  as: 

(1, 0, 0, 0, 0, 0, 0, 0)  means  Joe  is  broke  (Bill  has  seven  coins) 

(0,  1, 0, 0, 0, 0,  0, 0)  means  Joe  has  one  coin 

(0,0, 1,0, 0,0, 0,0)  means  Joe  has  two  coins 


(0,0, 0,0, 0,0, 0,1) 


means  Joe  has  seven  coins  (Bill  is  broke) 
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Here  the  "particle"  is  the  1  in  the  vector,  and  the  state  of  Joe's  wealth  at  the  close 
of  any  event  is  indicated  by  the  position  of  the  particle.  The  game  starts  with  the 
state  vector  being  V  =  (0,0,0, 1,0,0, 0,0),  i.e.,  with  Joe  having  three  coins,  and 
will  continue  through  a  sequence  of  vectors  indefinitely,  until  either  Joe  goes  broke 
(1 , 0, 0, 0, 0, 0, 0, 0)  or  Bill  goes  broke  (0, 0, 0, 0, 0, 0, 0, 1 ). 

If  the  game  is  just  about  to  start,  we  cannot  know  now  the  state  of  Joe's  wealth, 
say,  three  encounters  later.  However,  for  each  of  the  eight  states,  it  is  possible  now 
to  compute  the  probability  that  he  will  be  in  that  particular  state  at  the  close  of  the 
third  encounter.  To  make  this  and  other  computations,  first  a  transition  matrix  M  is 
displayed. 

1  0  0  0  0  0  0  0 

*0*00000 
0*0*0000 
00*0*000 

M  = 

000*0*00 
0  0  0  0  *  0  *0 

00000*0* 

0  0  0  0  0  0  0  1 

To  explain  the  construction  of  the  matrix,  it  is  noted  that  there  are  eight  rows 
and  eight  columns;  and  thus,  there  are  64  entries  within  the  matrix.  If  i  =  0,  1,  . . . , 
7,  and  j  =0,  1,  ...,  7,  and  p;j  are  the  entries,  then  any  p;:  appearing  in  the  i™ 
row  and  column  is  simply  the  answer  to  the  following  question:  "given  that  the 
particle  is  in  state,  wtat  is  the  probability  that  at  (immediately  after)  the  next 
encounter  (step)  the  particle  will  be  in  state  j?"  For  example,  if  at  some  instant 
Joe  has  six  coins,  then  as  an  immediate  result  of  the  next  matching  he  will  have 
either  five  or  seven  coins,  each  with  probability  1/2;  and  zero,  one,  two,  three, 
four,  or  six  are  impossible  (probability  is  zero).  These  facts  are  indicated  in  the 
next  to  last  row  of  M  by  the  ordered  entries  00000*0*. 

Now  to  show  the  use  of  M,  we  perform  the  multiplication  VQM  to  obtain  V] , 
the  state  vector,  just  after  the  first  matching. 

V,  =  VCM 


=(00*0*000) 

This  result  merely  shows  that  Jce  is  equally  likely  to  have  two  or  four  coins  after  the 
first  step,  and  it  is  impossible  to  have  at  that  time  any  other  sum. 
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Now,  performing  the  multiplication  Vj  M, 

V2  =  V,M 

V2  =  (o  I  o  i  0  1  0  0) 

which  simply  meons  that  at  the  end  of  the  second  step  Joe's  probabilities  are  as 
follows: 

P  (broke)  =  0 
P  (single  coin)  =  £ 

P  (two  coins)  -  0 
P  (three  coins)  =  £ 

P  (four  coins)  =  0 
P  (five  coins)  =  £ 

P  (six  coins)  =  0 
P  (seven  coins)  =  0 

Continuing, 

v3  =  v2m 

V3  =  (  }0  i  0  1  0  1  0) 

It  is  noted  that  at  each  step  the  sum  of  the  entries  in  the  state  vector  is  unity, 
as  required  by  the  fact  that  the  "particle"  must  at  all  times  be  somewhere. 

It  is  possible  to  continue  indefinitely  and  observe  that  as  the  number  of  steps 
increases  the  state  vector  will  approach  the  inevitable  result: 

=  (  f  0  0  0  0  0  0  f  ) 

This  simply  means  that  if  Joe  and  Bill  agree  to  play  until  on''  is  broke  the  probability 
that  Joe  goes  broke  eventually  is  . 
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A  MARKOVIAN  APPROXIMATION  TO  THIN-SLAB  ALBEDO 


Returning  now  to  the  problem  of  this  report,  it  is  noted  first  that  the  distance 
from  the  slab  face  to  a  neutron  impact  point  for  some  future  impact  of  interest  is  not 
exactly  amenable  to  the  Markov  approach  described  above,  because  this  distance  is 
a  continuous  variable  in  contradistinction  to  a  discrete  variable.  Thus,  the  number 
of  possible  states  is  uncountably  infinite.  Nevertheless,  it  is  proposed  to  obtain 
approximate  results  by  dividing  (conceptually)  the  slab  thickness  into  an  integer 
number  (^)  of  equally  thin  sub-slabs  (cells),  each  of  thickness  h.  This  gives  rise 
to  +  2  Markov  states  numbered  i  or  j  =  0,  1,  ^  +1. 

At  any  step  (serially  numbered  neutron  impact),  the  existence  of  the  neutron 
in  state  (or  cell)  zero  means  by  convention  that  it  has  already  returned  to  region 
Z  <  0  and,  thus,  will  never  actually  experience  that  numbered  impact.  At  any 
numbered  impact  (step),  state  five  means  that  the  neutron  has  experienced  that 
impact  at  some  point  such  that  4h  <  Z  <  5h  ^  T.  At  any  step  the  neutron  in  state 
(^+  1)  means  by  convention  that  it  has  already  passed  through  the  slab  of  thickness 
T  and  will  never  experience  that  or  later  impacts. 

The  reason  this  approach  will  not  be  exact  is  that  when  computing  a  typical 
element  m;  j  of  the  (-£+  2)  by  (  ^  +  2)  transition  matrix  M,  the  assumption  is  made 
that  if  the  last  impact  occurred  in  the  i^  cell,  it  occurred  at  a  midpoint  of  this 
cell.  This  tends  to  err  on  the  side  of  stability,  causing  the  diagonal  elements  mj| 
to  be  slightly  too  large.  It  is  expected,  however,  that  as  h  -  0,  i.e.,  as  the  number 
of  cells  is  increased,  the  computed  albedo,  as  well  as  other  desired  results,  will 
converge  rapidly  so  as  to  produce  results  with  sufficient  accuracy  for  practical 
purposes. 

If  by  convention  the  (-£+2)  elements  of  state  vector  Vj  are  the  probabilities 
that  the  first  impact  will  occur  in  the  corresponding  cell,  then  it  will  be  possible  to 
compute  Vi  exactly  (within  round-off  error). 

The  next  section  of  this  report  consists  of  a  specific  numerical  example.  There 
are  two  reasons  for  presenting  such  a  section.  First,  the  authors  find  that  it  is  easier 
to  explain  the  concepts  by  an  example  than  by  an  abstraction.  Second,  it  is  believed 
thnt  it  will  be  easier  to  follow. 


A  NUMERICAL  EXAMPLE 

In  order  to  achieve  a  better  understanding  of  the  physical  problem,  let  us 
consider  a  neutron  normally  incident  on  a  slab  whose  thickness  is  one  half  of  a  mean 
free  path.  Let  us  subdivide  the  slab  conceptually  into  five  layers  each  of  thickness 
h.  Assume  the  absorption  cross  section  is  one  fourth  as  large  as  the  scattering  cross 
section.  Then  we  have 
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®o  *  0 
h  ' 


i  =  0,  I,  5,  6 

The  "cells"  are  the  various  slab  regions.  For  example,  the  zeroth  cell,  or 
region,  is  the  void  which  the  incident  neutron  traverses  before  it  reaches  the  slab. 
The  first  cell  is  the  first  layer  of  the  slab  which  the  neutron  encounters  as  it  enters 
the  slab.  The  various  cells,  then,  are  defined  as  follows: 

i  =  0  Z  <  0 

i  =  1  0  <  Z  <  0.1 

i  =  2  0.1  <  Z  <  0.2 

i  =  3  0.2  <  Z  <  0.3 

i  =  4  0.3  <  Z  <  0.4 

i  =  5  ( I  =  5)  0.4  <  Z  <  0.5 

i  =  6  0.5  <  Z 

One  unit  of  distance  (for  example  along 
the  Z  axis)  is  equal  to  1/1^ . 

It  will  be  noted  that  a  point  on  a  cell  boundary  is  left  ambiguous,  and  the 
authors  without  apology  take  the  dodge  that  an  impact  will  occur  there  with  proba¬ 
bility  zero.  It  is  desired  to  approximate  the  fifteen  items  listed  in  this  report  under 
the  heading  "Statement  of  the  Problem." 

The  Vj  state  vector,  showing  the  position  by  cell  number  for  the  first  impact, 
will  now  be  computed.  The  probability  that  this  will  occur  in  cell  zero  is  zero, 
because  by  assumption  the  space  Z  <  0  is  vacuum  except  for  the  neutron  moving 
along  the  negative  Z  axis. 
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r 


■4 

_ i 


1 

10  Z, 


P  (Is*  impact  in  cell  1)  *  y 


-sec  0  I  Z 

sec0  e  0  *  dZ 


=  0.0951626 


! 


i 


5 

10  ^  -seceiZ 

P  (1st  impact  in  cell  5)  =  f  sec6o  Zf  e"^  °  *  dZ 

4 

10  z» 

=  0.0637894 


co 

.  n  -sec  6  L  Z 

P  (cell  6  entered  with  no  impact)  =  \  sec8e  I  e  * 

5 

'0l, 

=  0.6065307 

Thus,  the  transpose  of  state  vector  V|  is 


dZ 


vi  ■ 


0.0951626 

0.0861067 

0.0779125 

0.0704982 

0.0637894) 

0.6065307J 


This  simply  means  that,  a  priori,  the  position  of  the  first  impact  of  a  neutron 
w,th  a  nucleus  of  the  slab  mater, al  is  unknown,  but  nonetheless  the  probability  that 

Jj!'  ,OC7'  f°r  eX°mPre  ?n  ^  f°Urth  Cel1,  ?S  °-0704982  equivalently,  about 
7.05 /o  of  a  large  number  of  neutrons  will  experience  their  first  impact  somewhere 
within  cell  number  four. 
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As  a  prelude  to  the  construction  of  the  transition  motrix  M,  It  will  be  recalled 
that  for  an  arbitrary  directional  one-dimensional  W  coordinate  system  with  origin  at 
a  last  point  of  impact  from  which  an  unabsorbed  neutron  in  leaving  isotropically  (in 
three  space),  the  distribution  function  of  the  W  coordinate  of  the  next  impact  point 


Ei(ItW) 


W  <  0 


H(W)  = 


y  W  =  0 

-zfw 


Ei(-ZfW) 


W  >  0 


Let  us  now  measure  W  coordinates,  h,  and  T  in  distance  units  equal  to  (l/Z^.)  and  use 
the  above  expression,  a  table  of  the  exponential  function,  and  a  table  of  the  Ei(-) 
function  to  compute  Table  1  for  H(W)  versus  W. 

Table  1.  Tabulation  of  Values  of  Distribution  Function  H(W)  Versus  W 


W  H(W) 


W  H(W) 


W  H(W) 


W  H(W) 


-0.51 

-0.50 

-0.49 

-0.48 

-0.47 

-0.46 

-0.45 

-0.44 

-0.43 

-0.42 

-0.41 

-0.40 

-0.39 

-0.38 


0.1605531 
0.1633219 
0.1661514 
0. 1690435 
0.1719999 
0.1750229 
0.1781145 
0.1812770 
0.1845127 
0.1878240 
0.1912135 
0.1946840 
0.1982383 
0.2018794 


-0.37 

-0.36 

-0.35 

-0.34 

-0.33 

-0.32 

-0.31 

-0.30 

-0.29 

-0.28 

-0.27 

-0.26 

-0.25 

-0.24 


0.2056105 

0.2094350 

0.2133564 

0.2173784 

0.2215052 

0.2257046 

0.2300901 

0.2345576 

0.2391486 

0.2438687 

0.2487238 

0.2537203 

0.2588651 

0.2641657 


-0.23 

-0.22 

-0.21 

-0.20 

-0.19 

-0.18 

-0.17 

-0.16 

-0.15 

-0.14 

-0.13 

-0.12 

-0.11 

-0.30 


0.2696302 
0.2752676 
0.2810874 
0.2871003 
0.2933180 
0.2997535 
0.3064211 
0.3133370 
0.3205194 
0.327988 9 
0.3357693 
0.34388 77 
0.3523762 
0.3612725 


-0.09 

-0.08 

-0.07 

-0.06 

-0.05 

-0.04 

-0.03 

-0.02 

-0.01 

0.00 

+0.01 

+0.02 

+0.03 

+0.04 


0.3706221 
0.3804805 
0.3909176 
0.4020231 
0.4139172 
0.4267694 
0.4408360 
0.4565523 
0.4748353 
0.5000000 
0.525164 7 
0.5434477 
0.5591640 
0.5732306 
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At  each  impact  the  neutron  is  either  absorbed  permanently  then  and  there 
(with  probability  I  /I^  =  0.2),  or  it  moves  isotropically  in  three-space  from  there 
(with  probability  I^/I^  =  0.8).  To  assist  in  computing  the  mj  j  elements  of  M,  it  is 
convenient  to  orient  a  W  axis  parallel  to  the  Z  axis  with  W  increasing  with  increasing 
Z  and  with  the  origin  of  the  W  axis  at  the  center*  of  the  i^*  cell.  Then,  Table  1  is 
used  to  compute 


The  foregoing  seemingly  confusing  formulae  may  be  clarified  by  referring  to 
Figure  1,  which  is  a  geometric  display  of  the  cells  for  this  numerical  example.  The 
orientation  of  the  illustration  is  across  the  slab  perpendicular  to  the  faces  and  shows 
the  cell  locations  (Z  coordinates  of  the  edges  and  center  of  the  five  interior  cells). 
The  thickness  of  the  slab  is  0.5/Zj,. 


*  All  the  inaccuracy  of  this  method  stems  from  this  approximation;  i.e.,  the  mean 
effective  center  of  impact  for  the  cell  is  not  necessarily  at  the  center. 
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Cell  number 
Z  coordinate 


0  { 


1  !  2  ;  3  1  4  !  5  !  6 

-■ — r-^-T—r-n — r— r-r— i 

0.05  0.1  0.15  0.2  0.25  0.3  0.35  0.4  0.45  0.5 


T 


Figure  1.  Cell  number. 


Now,  for  example,  to  compute  m^j ,  referring  to  Figure  1,  it  is  noted  that  for 
the  neutron  to  move  from  cell  4  with  center  at  Z  =  0.35  to  cell  1,  it  is  necessary  to 
move  to  the  left  (decreasing  Z)  a  distance  someplace  between  0.25  and  0.35  units 
before  the  next  impact.  Now,  referring  to  Table  1,  !fi(-0.25)  -  H(-0.35)  =  0.2588651 
0.2133564  =  0.0455087.  But  for  any  impact  in  cell  4,  the  probability  of  immediate 
neutron  absorption  there  is  0.2 .  Thus, 

m..  =  0.8(0.0455087)  =  0.0364070 

4  I 

The  complete  transition  matrix  M  is  displayed  in  Table  2. 


Table  2. 


Transition  Matrix  M 


Z 

y-  =  0.8 


1.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.3311338 

0.1377324 

0.0747183 

0.0493234 

0.0364070 

0.0281935 

0.1424916 

0.2564155 

0.0747183 

0.1377324 

0.0747183 

0.0493234 

0.0364070 

0.1706851 

0.2070921 

0.0493234 

0.0747183 

0.1377324 

0.0747183 

0.0493234 

0.2070921 

0.1706851 

0.0364070 

0  0493234 

0.0747183 

0.1377324 

0  0747183 

0.2564155 

0.1424916 

0.0281935 

0.0364070 

0.0493234 

0.07 17183 

0.1377324 

0.3311338 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

1.0 

Using  V i  and  M  the  following  sequential  computations  are  performed: 


V]  AA 

V2M 

v3m 
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This  process  should  converge  reasonably  rapidly  to  Voo.al  imiting  vector  in  which 
the  zeroth  component  represents  total  albedo,  the  sixth  component  represents  proba¬ 
bility  of  transmission  through  the  slab,  and  the  other  components  are  zero.  For  any 
state  vector  V^,  one  minus  the  sum  of  the  components  is  the  probability  of  neutron 
absorption  before  the  k^1  impact.  Table  3  displays  the  results  of  the  first  few 
+  j  M  computations. 


Table  3. 


Successive  State  Vectors 


V] 

0.0 

0  0951626 

0.0861067 

0.0779125 

0.0704982 

0.0637894 

0  6065307 

v2 

0.0908481 

0.0277487 

0.0305911 

0.0302724 

0.0280093 

0.0237141 

0.6901224 

V3 

0.1223097 

0.0092891 

0.0107940 

0.0110863 

0.0104107 

0.0087482 

0.7206015 

V4 

0.1334728 

0.0032584 

0.0038411 

0.0040010 

0.0037865 

0.0031845 

0.7316297 

V5 

0.1374653 

0.0011608 

0.0013742 

0.0014388 

0.0013665 

0.0011506 

0.7356036 

v6 

0.1388972 

0.7370328 

Vao  (•*»•) 

0.13979 

0.0 

0.0 

0.0 

0.0 

0.0 

0.73792 

Estimated  probability  of  absorption  =  1  -  0.13979  -  0.73792  =  0.12229. 


In  the  process  of  computing  the  zeroth  component  of  any  j  vector,  using 
the  row  by  column  rule,  it  is  possible  to  note  in  detail  from  which  cells  the  increase 
in  total  albedo  springing  from  the  k"1  impact  came.  This  detail  is  recorded  in 
Table  4.  For  example,  the  probability  that  the  neutron  will  enter  the  zeroth  cell 
immediately  after  a  final  fourth  impact  in  cell  number  2  is  simply 

0.0038411  x  0.2564155  =  0.00098492 

where  0.003841 1  is  the  number  2  cell  component  of  V4  (Table  3),  and  0.2564155  is 
rr»20  read  from  Table  2.  The  result,  0.0009849,  is  recorded  in  Table  4.  The  last 
line  of  Table  4  contains  Q  values,  which  are  the  analysis  of  total  albedo  into  the 
cell  of  final  impact. 

In  a  similar  manner,  Table  5  is  a  record  of  entrances  into  cell  number  6, 
maintained  by  impact  number  and  cell  source. 

Total  albedo  (probability  of  entering  cell  zero),  or  proportion  of  a  very  large 
number  of  incident  neutrons  which  backscatter,  consists  of  those  which  have  exper¬ 
ienced  exactly  one  impact,  exactly  two  impacts,  exactly  three  impacts,  etc.  From 
the  left  column  of  Table  3,  by  taking  differences,  it  is  possible  to  construct  Table  6, 
which  is  a  breakdown  of  neutrons  entering  cell  zero  into  the  number  of  impacts 
before  returning  to  the  region  Z  <  0.  Of  course,  it  is  computationally  expedient 
during  the  matrix  multiplication  leading  to  Table  3  to  record  entries  in  Table  6. 

(This  avoids  the  need  for  actual  subtraction.) 
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Table  4.  Probability  of  Neutron's  Entering  Zeroth  Cell,  Having 
Had  i  Impacts  With  i™  Impact  Occurring  in  Cell  j 
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0.6065307  0.01960  0.02280  0.02601  0.02945  0.03354  0.73793 


Table  6. 


Number  of  Impacts  Analysis  of 
Albedo  (Zeroth  Cell  Entrants) 


0.8 


Number  of  Impacts 

Probability 

1 

0.0908481 

2 

0.0314616 

3 

0.0111631 

4 

0.0039925 

5 

0.0014319 

Total  of  5 

0.1388972 

6  and  more  (est.) 

0.00089 

Estimated  total  albedo 

0.13979 

In  a  similar  manner,  the  estimated  number  of  historic  impacts  for  neutrons 
entering  cell  6  (passing  through  the  slab  never  to  return)  are  listed  in  Table  7. 


Table  7. 


Number  of  Impacts  Analysis  of 
Neutrons  Transmitted  Through  Slab 


Number  of  Impacts 

Probability 

0 

0.6065307 

1 

0.0835917 

2 

0.0304791 

3 

0.0110282 

4 

0.0039739 

5 

0.0014293 

Total  0  to  5  inclusive 

0.7370328 

6  and  more  (est.) 

0.00089 

Estimated  total  through  slab 

0.73792 
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In  order  to  break  down  the  0.12229  absorbed  (beneath  Table  3)  into  the  impact 
number  (i.e.,  first,  second,  third,  etc.)  at  which  absorption  takes  place,  Table  8  is 
constructed  by  first  summing  the  components  of  the  successive  state  vectors  (Table  3) 
and  then  subtracting  to  find  the  amounts  of  absorption  at  successive  impacts. 


Table  8.  Absorption  of  Neutrons  by  Cell  and  Impact  Number 


Impact 

Number 

(i) 

6 

2  V,, 

j  =  0 

Absorbed 

Cell  1 

Cell  2 

B 

1 

Cell  4 

Cell  5 

1 

1.0000000 

0.0786939 

0.0190325 

0.0172213 

0.0155825 

0.0127579 

2 

0.9213061 

0.0280666 

0.0055497 

0.0061  B2 

0.0047428 

3 

0.8932395 

0.0100657 

0.0018578 

0.0021588 

0.0022173 

0.0020821 

0.0017496 

4 

0.8831739 

0.0036143 

0.0006517 

0.0007682 

5 

Est  6,  7,  8,  .  .  . 
Est  ao 

0.8795596 

0.87771 

0.0012981 

0.00055 

0.0002322 

0.00011 

0.0002748 

0.00011 

0.0002733 

0.00011 

0.0002301 

0.00011 

Total 

0.12229 

0.02743 

0.02665 

0.02505 

0.02292 

0.02023 

The  foregoing  computations  have  been  concerned  in  the  main  with  proportion 
absorbed,  proportion  passing  through,  and  total  albedo.  It  is  now  time  to  come  to 
grips  with  the  problem  of  analyzing  the  total  albedo  with  reference  to  the  directional 
(0)  density,  g(0),  of  backscatter  and,  of  even  more  interest  to  the  physicist,  the 
differential  albedo  a(0,  0)  for  0  <  0  <■  Tr/2 .  It  is  essential  to  note  that  these  differ 
from  each  other. 

First,  define  Cl  ,  the  total  albedo,  to  be  the  probability  that  eventually  the 
neutron  backscatters  to  Z  <  0 .  This  is,  of  course,  a  function  of  L ,  T,  and  0  . 

a  S f  r  O 

Then  ol,  the  differential  albedo  per  steradian  is  defined  so  that 

\  a  d  n  =  U 
BHS 

where  BHS  is  the  back  hemisphere  in  the  region  Z  <  0,  with  base  in  the  plane  Z  =  0, 
and  d  Q  is  the  differential  solid  angle  so  that 

d  0  =  sin  0  d0  d<P 
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The  defining  relation  can  be  rewritten  as 


2v  ir/2 

n  a  sin  Q  66  d<p  =  (? 

0  0 

In  general,  a  is  a  function  of  both  6  and  4> .  For  the  isotropic  case  of  this 
report,  a  is  independent  of  <p ;  nevertheless,  it  is  well  to  remember  that  a/ &  is  a 
two-dimensional  directional  density  of  the  backscattered  neutron.  After  the  first 
integration  is  performed, 

f 

2n  j  asin  9  d0  =  (i 

0 


or 


de  =  a 


where 


2 

g(9)  =  -jj-  tr  asin  0 

a  density  which  displays  the  relative  frequency  of  return  at  all  the  values  of  0, 

o  <  0  <  it/: 2 . 

For  a  neutron  that  is  backscattered,  having  had  its  last  impact  at  a  known 
depth  Z,  the  angles  (0,  <p)  of  return  are,  of  course,  random  variables  with  some  joint 
density,  say  f (<p,  8),  where 


\  ffo,  0)dO  =  1 
BHS 

If  an  impact  is  at  a  depth  Z  and  the  direction  of  a  nonabsorbed  neutron  is 
0  £  0  <  it/ 2 ,  0  £  <p  <  2 IT,  then  the  probability  that  it  will  leave  the  slab  (via  Z  =  0) 
with  no  more  impacts  is  e“^t  Zsec0  ancj  because  of  isotopy,  f (<f>,  0)  must 

be  proportional  to  this.  Restated  with  proportionality  constant  K, 
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-I  2  sec  0 

f(0,0)dO  =  Ke  T  dO 


To  evaluate  K, 


-rr- 

o  o 


-I  Z  sec  0 

o  sin  6  d0d0 


‘  I  •  ’ 


-I.  Zsec  6 


sin  0d0 


T“ztz  1 

=  2tr  K  e  +  Z  Ei(-If  Z) 


4trH(-ItZ) 


where  H  is  defined  as 


i  r_itz  i 

(-Z,  Z)  =  y  e  +  I,  Z  Ei(-I,  Z) 


The  density  function  then  becomes 
-Ij  Z  sec  0 

f(0,0)  =  4ffH(-It  Z) 

This  is  the  directional  density  of  return  for  a  neutron  returning  with  its  final 
impact  having  been  at  depth  Z.  The  bottom  of  Table  4  estimates  the  breakdown  of 
total  albedo  into  the  cell  of  final  impact  and  thus  gives  depths  Z  roughly  (within  a 
half  cell  width).  If  it  is  now  assumed  that  a  neutron  leaving  a  cell  had  its  impact 
at  the  middle  of  the  cell,  then  a  rough  estimate  of  a  may  be  had  by  summing  over 
all  cells.  Ideally,  if  <P( Z)  dZ  is  the  a  priori  probability  that  an  incident  neutron 
will  eventually  backscatter  from  a  last  impact  at  depth  Z,  then 


a  (9,0)  =  j*  0(Z)f(0,  6)dZ 
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This  will  now  be  approximated  roughly.  The  approximation  moy  be  improved 
ultimately  by  dividing  the  slab  into  more  cells.  The  approximation  will  be  made  in 
Table  9  for  various  angles  0.  To  do  this,  f(0,  0)  is  computed  and  entered  into  the 
third  column  for  various  values  of  0 .  The  Q  values  of  the  fourth  column  have  been 
taken  from  the  last  line  of  Table  4.  The  sum  of  the  products,  Qj  f (<t>,  0),  then  is  an 
estimate  of  differential  albedo.  The  results  in  Table  9  are  plotted  in  Figure  2. 

Using  the  results  of  Table  9,  the  density  g(0)  is  computed  in  Table  10  and 
plotted  in  Figure  3  as  (J  g(0).  The  area  under  this  curve  is  0.1398,  the  estimated 
total  albedo.  It  is  noted  that  this  function  has  a  mode  at  about  0  =  58  degrees. 

This  simply  means  that  a  greater  percentage  of  the  total  albedo  scatters  back  into 
a  neighborhood  d0  about  an  angle  0  of  approximately  58  degrees  than  into  other 
neighborhoods  of  the  same  size.  This  does  not  say  that  differential  albedo  is 
maximum  near  58  degrees;  differential  albedo  has  its  maximum  at  0  =  0  (Figure  2). 

Computations  for  the  scatter  through  angular  distribution  similar  to  those 
for  albedo  (Tables  9  and  10  and  Figures  2  and  3)  could  be  made  using  the  same 
technique  —  a  novel  feature  being  that  here  fl/2  <0  i  if.  These  computations  were 
not  actually  performed  for  this  example.  The  density  function  for  total  neutrons 
going  through  would  contain  the  additional  novel  feature  that  there  is  a  spike  of 
mass  (in  this  case  0.60653/0.73792)  in  the  direction  -0o  =  TT,  corresponding  to  the 
probability  of  passing  directly  through  the  slab  with  zero  impacts. 

There  remains  the  task  of  estimating  the  expected  number  of  neutron  impacts 
and  to  analyze  this  into  various  classifications.  There  are  two  methods  (say,  A  and 
B)  available  for  estimating  the  expected  number  of  impacts. 

To  explain  method  A,  consider  Table  3.  For  Vj  the  sum  of  the  components 
1  through  5  (0.3934693)  is  the  probability  of  experiencing  a  first  impact.  The 
probability  of  experiencing  a  second  impact  is  the  sum  of  the  components  1  through  5 
ofN^.  Thus,  if  these  component  sums  are  summed  over  all  Vj ,  ?  =  1,  2,  3,  ...,  this 
double  summation  will  be  an  estimate  of  the  expected  number  of  total  impacts. 

Method  B  is  a  direct  application  of  the  definition  of  expected  value,  namely 

ao 

hP  (exactly  h  impacts)  =  E  (number  of  impacts) 
h=1 


To  use  this  formula  it  will  be  necessary  to  estimate  P  (exactly  h  impacts)  as  the  sum 
of  the  following  three  quantities:  (1)  probability  of  entering  cell  0,  having  had 
exactly  h  impacts;  (2)  probability  of  entering  cell  6,  having  had  exactly  h  impacts; 
and  (3)  probability  of  being  absorbed  at  the  h*^1  impact. 

By  the  nature  of  the  convergences,  method  A  should,  in  practice,  be  the  more 
accurate,  and  method  B  will  constitute  an  approximate  check. 
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Table  9.  Computation  of  Differential  Albedo 


5tt/18 


*.z 

f(0,  6) 

0 

o  ffo,  e) 

mm 

0.1828782 

0.04544 

0.0083100 

0.2136937 

0.03421 

0.0073105 

■H' 

0.2394104 

0.02598 

0.0062199 

0.2628340 

0.01976 

0.0051936 

0.45 

0.2848776 

0.01442 

0.0041079 

0.1822923 

0.2116464 

0.2355999 

0.2569960 

0.2767680 


a (9/  0)  = 


0.04544 

0.03421 

0.02598 

0.01976 

0.01442 


a  (<f>,  ff/9)  = 


n 

0.1801068 

0.04544 

0.2041250 

0.03421 

■ss 

0.2218114 

0.02598 

0.2361884 

0.01976 

0.45 

0.2482972 

0.01442 

0.1778666 

0.1966027 

0.2083561 

0.2163763 

0.2218461 


a(<p,  ir/9)  = 


0.04544 

0.03421 

0.02598 

0.01976 

0.01442 


0.0311419 


0.0082834 

0.0072404 

0.0061209 

0.0050782 

0.0039910 


0.0307139 


0.0081841 

0.0069831 

0.0057627 

0.0046671 

0.0035804 


0.0291774 


0.0080823 

0.0067258 

0.0054131 

0.0042756 

0.0031990 


a  (6,  5it/18)  =  0.0276958 


Continued 


Table  9.  Continued 


6 

r.2 

H4>,  S) 

O 

Q  f (6,  0) 

it/ 3 

■S9 

0.1739591 

0.04544 

0.0079047 

0.1839279 

0.03421 

0.0062922 

0.1864530 

0.02598 

0.0048440 

■19-: 

0.1852160 

0.01976 

0.0036599 

0.45 

0.1816460 

0.01442 

0.0026193 

a  (6,  ff/3)  = 

0.0253201 

7ir/ 18 

KS 

0.1661066 

0.04544 

0.0075479 

0.1601279 

0.03421 

0.0054780 

0.1480024 

0.02598 

0.0038451 

0.1340471 

0.019/6 

0.0026488 

0.45 

0.1198628 

0.01442 

0.0017284 

a  (6/  7ff/18)  = 

0.0212482 

4n/9 

0.05 

0.1441539 

0.04544 

0.00655 04 

0.15 

0.1046612 

0.03421 

0.0035805 

0.25 

0.0728561 

0.02598 

0.0018928 

0.35 

0.0496 974 

0.01976 

0.0009820 

0.45 

0.0334688 

0.01442 

0.0004826 

a  (6,  4tt/9)  = 

0.0134883 

1 7rr/36 

1 

0.1083248 

0.04544 

0.0049223 

sim 

0.0444110 

0.03421 

0.0015193 

0.25 

0.0174572 

0.02598 

0.0004535 

0.35 

0.0067243 

0.01976 

0.0001329 

0.45 

0.0025571 

0.01442 

0.000369 

a  (6,  1 7tr/36)  = 

0.0070649 
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Table  10.  The  Density  g(0), 


0 

2trsin  0 

a 

fig®) 

g(0) 

0.0 

0.0 

0.0311419 

0.0 

0.0 

ff/9 

2.148976 

0.0307139 

0.06600 

0.47214 

2  ir/9 

4.038754 

0.0291774 

0.11784 

0.84298 

5ir/ 18 

4.813199 

0.0276958 

0.13331 

0.95364 

it/ 3 

5.441398 

0.0253201 

0.13778 

0.98562 

7it/18 

5.904263 

0.0212482 

0.12545 

0.89742 

4tt/9 

6.187730 

0.0134883 

0.08346 

0.59704 

17tr/36 

6.259276 

0.0070649 

0.04422 

0.31633 

ir/2 

6.283185 

0.0 

0.0 

0.0 

Proceeding  with  method  A,  Table  11  is  prepared  from  Table  3  by  computing 
the  sums 


5 


j  =  l 


and  then  finally  estimating  the  sum  of  all  these,  i=  1,  2,  . . .,  obtaining  the  result 
0.612  for  the  expected  total  number  of  impacts  in  the  offing  for  a  neutron  directed 
at  the  slab. 

Proceeding  with  method  B,  Table  12  is  constructed  by  transcribing  Tables  6, 

7,  and  8  to  columns  2,  3  and  4  of  Table  12,  summing  these  to  form  the  column 
"Probability  of  exactly  h  impacts,"  and  then  forming  the  last  column  by  the  indicated 
multiplication. 

The  reason  the  result  of  Table  12  is  intrinsicly  less  accurate  than  that  of 
Table  11  is  that  the  figure  0.0150  in  the  lower  right  corner  of  Table  12  is  a  rough 
infinite-sum  estimate,  made  rough  because  the  ratios  of  successive  entries  in  the 
last  column  stabilize  much  less  rapidly  thjn  those  of  Table  1 1. 

A  breakdown  of  the  expected  0.61233  impacts  can  be  approximated  as  shown 
in  Table  13.  The  expected  number  of  absorbing  impacts  is  computed  simply  as 
0,2  x  0.61233  -  0.12247.  To  compute  the  expected  number  of  impacts  invested  in 
absorbed  neutrons,  the  cross-products  of  the  first  and  third  column  of  Table  8  are 
summed.  Here  again  the  slow  convergence  introduces  inaccuracy  in  the  estimate 
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for  i  =  6,  7,  ....  Similarly,  Tables  6  and  7  are  used  to  estimate  the  expected 
number  of  impacts  of  the  incident  neutron  to  be  invested,  respectively,  in  albedo 
and  scatter-through.  The  Table  13  round-off  error,  0.61310  versus  0.61233,  could 
be  reduced  by  extending  Table  3  with  some  more  matrix  multiplications,  computing, 

saY/  v6'  V7,  and  VR- 

If  a  neutron  Is  backscattered,  the  expected  number  of  impacts  it  has  is 

E(|sj  )  =  -  =  ]  544 

'  b;  0.13979  ^ 

If  a  neutron  has  been  absorbed,  the  expected  number  of  impacts  it  has  had  is 

E<V  -  01222?  -  ,'558 


If  a  neutron  is  scattered  through,  the  expected  number  of  impacts  it  has  had  is 

r,M  x  0.20641 

E(Ns>  -  0J3139  -  ’•571 


Table  11. 


Expected  Number  of  Impacts 


* 


• 

1 

impact) 

1 

0.3934693 

2 

0.1403356 

3 

0.0503283 

4 

0.0180714 

5 

0.0064908 

6 

0.002331  (est.) 

7 

0.000837  (est.) 

8,  9,  . . . 

0.00047  (est.) 

00 

2  • 

0.61233  (est.) 

<  =  > 

*  Probability  of  experiencing  the  i^1  impact  Is  computed  as  2  v;j  (from  Table  3). 

j  =  1 

Estimated  expected  number  of  impacts  is  0.61233. 
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Table  12. 


Check  Estimate  of  Expected  Number  of  Impacts 


h 

0) 

With  Exactly  h  Impacts 

(2) 

Probability 
of  Exactly 
h  Impacts 

(1)  •  (2) 

P(Backscatter) 

P(Scatter- 

Through) 

P(A’  orbed) 

1 

0.0908481 

0.CS35917 

0.0786939 

0.2531337 

0.2531337 

2 

0.0314616 

0.0304791 

0.0280666 

0.0900073 

0.1800146 

3 

0.0111631 

0.0110282 

0.0100657 

0.0322569 

0.0967707 

4 

0.0039925 

0.0039739 

0.0036143 

0.0115807 

0.0463227 

5 

0.0014319 

0.0014293 

0.0012981 

0.0041593 

0.0207965 

6,  7,  ...  Est 

0.00089 

0.00089 

0.00055 

0.00233 

0.0150 

Expected  total  number  of  impacts  0.6120  Est 


Table  13. 


Analysis  of  Expected  Total  Number  of  Impacts 


Expected  number  of  absorbing  impacts  0.12247 

Expected  number  of  impacts  of  neutron  before  absorbing  impact  0.06810 


Expected  number  of  impacts  invested  in  absorption  0.19057 

Expected  number  of  impacts  invested  in  albedo  0.21612 

Expected  number  of  impacts  invested  in  scatter-through  0.20641 


Total  expected  impacts  0.61310 


THE  PERFECT  SCATTERING  MEDIUM 

For  the  numerical  example  solved  in  the  preceding  section,  it  will  be  recalled 
that  =  0.8;  i.e.,  the  slab  material  could  absorb  neutrons,  as  well  as  scatter 

them.  The  next  case  to  be  considered  is  the  case  of  the  perfect  scatterer.  That  is, 

*5  =  *t 
Za  =  0 
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As  before,  we  choose 


h  '  10  Zf 

®o  =  0 


j  .  0,  1, 2, 3, 4, 5, 6 

Since  the  computations  are  similar,  only  results  will  be  presented.  Tables  2A 
through  7 A,  corresponding  to  Tables  2  through  7,  give  results  for  the  perfect  scattering 
case.  Figure  2A  displays  differential  albedo  for  the  perfect  scattering  case. 

The  expected  number  of  collisions  of  a  neutron  before  it  leaves  the  slab  is 
found  to  be 

E(N)  =  0.711 


The  expected  number  of  collisions  of  a  neutron  which  has  backscattered  from 
the  slab  is 

E(Mb)  =  1.792 

The  expected  number  of  collisions  of  a  neutron  which  is  scattered  through  the 
slab,  i.e.,  has  suffered  at  least  one  impact  while  penetrating,  is 

E(N$)  =  1.824 

The  expected  number  of  collisions  of  a  neutron  which  penetrates  the  slab,  with 
or  without  collisions,  is 

E(Nf)  =  0.437 


Table  2A. 


Transition  Matrix 


l 

0.4139172 
0.3205194 
0.2588651 
0.2133564 
0.1781 145 
.0 


0 

0.1721655 

0.0933979 

0.0616543 

0.0455087 

0.0352418 

0 


0 

0.0933979 

01721655 

0.0933979 

0.0616543 

0.0455087 

0 


0 

0.0616543 

0.0933979 

0.1721655 

0.0933979 

0.0616543 

0 


0 

0.0455087 

0.0616543 

0.0933979 

0.1721655 

0.0933979 

0 


0 

0.0352418 

0.0455087 

0.0616543 

0.0933979 

0.1721655 

0 


0 

0.1781145 

0.2133564 

0.2588651 

0.32C5194 

0.4139172 

1 
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Table  3A.  Successive  State  Vectors  (Z  =  I ) 

s  t 


Vi 

0 

0.0951626 

0.0861066 

0.0779125 

0.0704982 

0.0637894 

0.6065307 

V2 

0.1135602 

0.0346859 

0.0382389 

0.0378405 

0.0350116 

0.0296427 

0.7110203 

V3 

0.1627189 

0.0145142 

0.0168657 

0.0173224 

0.0162667 

0.0136691 

0.7586440 

V4 

0.1845217 

0.0063641 

0.0075021 

0.0078144 

0.0073955 

0.0062196 

0.7801833 

V5 

0. 1942691 

0.0028339 

0.0033548 

0.0035126 

0.0033362 

0.0028090 

0.7898852 

0.1986388 

0.7942470 

voo  (•»*■) 

0.2021809 

0 

0 

0 

0 

0 

0.7978191 

Table  4A.  Probability  of  Neutron's  Entering  Zeroth  Cell  Having  Had 
i  Impacts  With  i^1  Impact  Occurring  in  Cell  j  (I  =  Z^) 


j 

Total 

BIS 

1 

2 

3 

4 

5 

0.0393894 

0.0275988 

0.0201688 

0.0150412 

0.0113618 

0.0143571 

0.0122563 

0.0097956 

0.0074699 

0.0052798 

0.0060077 

0.0054053 

0.0044842 

0.0034706 

0.0024347 

0.0026342 

0. 0024046 

0.0020229 

0.0015779 

0.0011078 

5 

0.0011730 

0.0010753 

0.0009093 

0.0007118 

0.0005003 

5 

-  Pij 
i  =  l 

1 - 

0.0635614 

0.0487408 

0.0373S08 

0.0282714 

0.0206844 

0.1986388 

Est  >.  Pij  =  Qj 
i  =  1 

0.0645274 

0.0496263 

0.0381053 

0.0288349 

0.0210869 

0.2021809 
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Table  5A.  Probability  of  Neutron's  Entering  Sixth  Cell,  Having  Had 
i  Impacts  With  i^1  Impact  Occurring  in  Cell  j  (1^  =  1^) 


■ 

j 

T  otal 

■ 

0 

1 

2 

3 

4 

5 

n 

0.6065307 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0169498 

0.0183714 

0.0201688 

0.0225960 

0.0264035 

0.0 

0.0061780 

0.0081585 

0.0097956 

0.0112219 

0.0122696 

0.0 

0.0025852 

0.0035984 

0.0044842 

0.0052138 

0.0056579 

0.0 

0.0011335 

0.0016006 

0.0020229 

0.0023704 

0.0025744 

L5 

0.0 

0.0005048 

0.0007159 

0.0009093 

0.0010693 

0.0011627 

5 

iV” 

0.6065307 

0.0273514 

0.0324447 

0.0373808 

0.0424714 

0.0480681 

0.7942470 

OO 

Est  X  Pij 

i  =6 

0.0005205 

0.0006174 

0.00071 13 

0.0008082 

0.0009147 

0.0035721 

oo 

Est  X  Pij 

i  =  J 

0.6065307 

0.0278719 

0.0330621 

0.0380921 

0.0432796 

0.0489828 

0.7978191 

Table  6A.  Number  of  Impacts 
Analysis  of  Albedo 
(Zeroth  Cell  Entrants) 


Number  of  Impacts 

Probabi  1  ity 

1 

0.1135602 

2 

0.0491587 

3 

0.0218029 

4 

0.0097473 

5 

0.0043697 

Total  of  5 

0.1986388 

6  and  mors  (est.) 

0.003542 

Estimated  total  albedo 

0.202181 

Table  7A.  Number  of  Impacts 
Analysis  of  Neutrons 
Transmitted  Through 
Slab  (I  =  If) 


Number  of  Impacts 

Probobi  lity 

0 

0.6065307 

1 

0.1044896 

2 

0.0476236 

3 

0.0215394 

4 

0.0097019 

5 

0  0043618 

Total  of  5 

0.7942470 

6  and  more  (est.) 

0.003572 

Estimated  total  transmission 

0.797819 
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scattering  medium. 


SUMMARY 


The  two  cases  that  have  been  considered  have  both  involved  a  neutron  normally 
incident  on  a  plane  slab  with  a  thickness  of  one-half  mean  free  path.  Using  a  discrete- 
ordinate  approximation,  the  neutron-diffusion  problem  has  been  solved  by  employing 
a  Markov  chain  process  to  represent  the  random  walk  phenomenon  of  isotropic  scatter¬ 
ing,  with  or  without  absorption. 

The  15  questions  posed  in  the  section  Statement  of  the  Problem  have  been 
solved,  as  reviewed  below: 


1.  For  Z  /Z.  =  0.8,  the  total  albedo  (3  was  estimated  to  be  0.13979.  For 
I  =  I. ,  (3  =  0.20218. 

2.  The  estimated  differential  albedo  a  for  the  two  cases  is  shown  in  Figures  2 
and  2A. 

3.  The  analyses  of  total  albedo  into  number  of  impacts  are  presented  in 
Tables  6  and  6A. 

4.  The  analyses  of  total  albedo  into  depth  probability  density  of  last  impact 
before  backscatter  from  the  slab  are  given  in  Tables  4  and  4A. 

5.  The  probability  of  absorption  for  the  case  2^/Zj.  =  0.8  was  estimated  to  be 
(3  =  0.12229.  For  the  perfect  scattering  case,  clearly  the  probability  of  absorption 


is  zero. 


6.  The  analysis  of  (3°  into  number  of  impacts  for  the  case  where  Zj/Ij,  =  0.8 
is  shown  in  Table  8. 

7.  The  analysis  of  (3°  into  depth  density  of  absorption  point  for  the  case  where 
Is/Ij.  =  0.8  is  also  shown  in  Table  8. 

8.  The  probabilities  of  transmission,  with  or  without  collision,  were  estimated 
to  be  0.7392  for  Z  /Zj,  =  0.8  and  0.79782  for  Z$  =  Zj, . 

9.  The  analyses  of  transmitted  neutrons  into  number  of  collisions  are  given  in 
Tables  7  and  7A. 

10.  The  analyses  of  neutrons  scattered  through  into  depth  density  of  last 
impact  before  penetration  are  given  in  Tables  5  and  5A. 

11-15.  The  expected  numbers  of  impacts  were  estimated  to  be 


w 

0.8 

^*t 

E(N) 

0.612 

0.711 

E(Nb) 

1.546 

1.792 

E(NS) 

1.571 

1.824 

E(N  ) 

0.280 

0.437 

E(N0) 

1.558 

— 

where  the  subscripts  b,  s,  t,  and  o  refer  to  backscatter,  scatter  through,  transmission, 
and  absorption,  respectively. 
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All  of  the  numerical  values  obtained  are  certain  to  be  somewhat  in  error 
because  of  the  approximation  of  assuming  that  all  collisions  occur  at  midpoints  of 
the  slab  layers.  The  approximation  can  be  improved  by  increasing  the  number  of 
subdivisions  or  layers.  A  consequence  will  be  that  the  transition  matrix  will  need 
to  become  larger. 


PUNS  FOR  THE  FUTURE 

It  is  planned  to  devise  and  employ  a  computer  program  to  perform  calculations 
with  increasingly  fine  slab  layers.  Results  will  be  examined  for  convergence.  If 
results  for  infinitely  thin  layers  can  be  approximated  by  extrapolation,  the  problem 
will  be  solved  to  a  high  order  of  approximation. 

The  computer  program  will  also  be  used  to  generate  albedo  data  for  slabs  of 
various  thicknesses,  for  different  values  of  Z  /Z^,  and  for  various  polar  angles  of 
neutron  incidence. 

Finally,  attempts  will  be  made  to  treat  more  difficult  problems,  such  as  the 
case  where  neutron  scattering  is  not  isotropic,  and,  if  possible,  the  case  where  Ij 
and  are  functions  of  the  neutron  energy,  which  changes  with  each  scattering 
event. 

Albedo  values  calculated  by  the  Markov  matrix  method  will  be  compared  with 
available  experimental  values.  If  agreement  is  satisfactory  to.  known  cases,  albedo 
values  can  be  generated  for  use  in  shielding  problems  where  backscatter  is  important. 
Also,  the  method  can  readily  be  generalized  to  treat  cases  of  penetration  through 
thin  slabs  of  scattering  and  absorbing  shields. 
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GLOSSARY  OF  SYMBOLS 


Total  albedo;  i.o.,  probability  that  ovontually  tho 
nautron  backscotters  out  of  tho  slob.  This  is,  of 
course,  a  function  of  2„,  21(  T,  and  0O  , 

Total  absorption;  i.o.,  probability  that  ovontually 
tho  nautron  is  absorbod.  This  is,  of  course,  a 
function  of  2g,  2t,  T,  and  0O  . 

Total  transmission;  i.o.,  probability  that  ovontually 
tho  nautron  is  transmitted  through  tho  slab.  This  is, 
of  course,  a  function  of  2^,  21,,  T,  ond  0O  . 

Probability  that  tho  neutron  passes  through  tho  slab 
with  no  impacts;  i.o.,  uncollided  transmission. 

=  .*T^soc«o. 

Back  hemisphere  in  region  2  <  0,  with  base  in 
plane  Z  =  0  . 

Differential  solid  angle,  dfi  =  sin  8  d8d<£. 


N0  If  the  neutron  is  absorbed,  this  represents  the  number 
of  impacts  it  has  experienced.  Thus,  Na  >  )  . 

Ns  If  evntually  the  neutron  scatters  through  the  slab, 
this  represents  the  number  of  impacts  it  hos  exper¬ 
ienced.  Thus,  N,  >  1  . 

N(  If  eventually  the  neutron  goes  through  the  slob,  this 

represents  the  number  of  impacts  it  has  experienced. 
Thus,  Nt  >  0. 

Pij  Probability  of  entering  zeroth  (sixth)  cell,  having  had 
i  impacts  with  the  i*h  impact  occurring  in  cell  j. 

P(-)  Probability  of  (■)■  For  example,  P(j)  denotes  proba¬ 
bility  that  the  neutron  will  experience  exoctly  j 
impacts. 

Qi  Probability  of  backscattering  with  the  last  impact 

in  the  jth  cell. 


Expected  value  operator.  Thus,  E(N)  simply  denotes 
the  expected  (or  average)  number  of  impacts  to  be 
experienced  by  the  neutron,  where  average  is  taken 
over  an  infinite  number  of  trials  of  the  experiment: 

"shoot  the  neutron  at  the  slab  under  the  given 
conditions." 

Directional  (<f>,  8),  two-dimensional,  density  of  the 
backscotter  from  a  fixed  known  depth,  Z,  where 

Directional  (0),  one-dimensional,  density  of  back- 

scatter,  where  f17  ^  g(8)d8  =  1. 

0 

Sub-slab  thickness;  so  that  T/h  is  a  positive  integer 

number  of  sub-slobs,  all  of  equal  thickness;  and  the 

number  of  Markov  states  is  (2  +  T/h).  The  symbol  h  X, 

is  also  used  as  an  index  for  counting  impacts. 

Probability  density  function  of  W  projection  of  next 
impact. 

Probability  distribution  function,  where  dH(W)/dW  = 
h(W) . 

Markov  state  indices,  where  i,  j  =  0,  1,  2,  ...  ,  T/h , 

T/h  +  1  .  Also,  i,j  ore  used  to  count  impacts. 

A  proportionality  constant  used  to  compute  f(d>,  8)  . 

Entry  in  M  ot  i*^  row  and  j**1  column. 

Transition  matrix. 

Total  number  of  neutron  impacts.  Thus,  N  ^  0  . 

If  eventually  the  neutron  bockscotters  to  Z  <  0, 
this  represents  the  number  of  impacts  it  has 
experienced.  Thus,  Nfc  )  . 


S  Distonce  (scaler)  measured  in  same  units  as  the 

Cartesian  coordinates. 

T  Known  slab  thickness.  Thus,  space  occupied  by  the 

slab  is  (-00  <  X  <  oo),  (-CD  <  Y  <  oo),  (0  <  Z  T) . 

Vjj  j*h  cell  entry  of  i*h  state  vector. 

V  2f  W  sec  /3  . 

Vj  i*h  state  vector. 

W  An  arbitrarily  directed  one-dimensional  coordinate 

system  in  the  same  distonce  units  as  the  Cartes  ion 
system,  and  with  origin  of  on  impact  point  of  the 
neutron. 

,  Z  A  fixed  set  of  right-handed  Cartesian  axes  with  the 
initial  slob  foce  being  the  (X,  Y)  plane  and  the  -Z 
axis  being  totally  outside  the  slab. 

<2.  Differential  albedo  per  steradian;  a  function  of 

^lo,  2,,  T,  3j,  and  8,  where  the  integral  over  the 

backward  (Z  <  0)  hemisphere  /  adf 2  =  (?  or, 

•  2rr  *2  BM5  & 

equivalently,  Jg  Jg  a  sin  8d$  d<f>  =  (I  . 

a!  Differential  scatter-through  (collided  transmission) 

per  steradian;  o  function  of  2^,,  2,,  T,  90,  ond  8, 
where  the  integral  over  the  forward  (Z  >  T)  hemi¬ 
sphere  jpyis  a,dfl  =  (?-(?„  or,  equivalently, 

f2n  fn  a,nin8d8d4>  =  (?  -  3o  • 

J0 

/S  Angle  between  a  roy  and  the  positive  W  axis. 

7)  Angle  between  projection  of  ray  onto  a  plane  perpen¬ 

dicular  to  the  W  axis  ot  W  =  0  and  a  fixed  direction 
in  this  plane. 
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If  th*  neutron  i*  absorbed,  this  represents  the  number 
>f  impacts  it  has  experienced.  Thus,  N0  >_  1  . 

If  eventually  the  neutron  scatters  through  the  slab, 
this  represents  the  number  of  impacts  it  has  exper¬ 
ienced.  Thus,  Nt  >  1  . 

if  eventually  the  neutron  goes  through  the  slab,  this 
'epresents  the  number  of  impacts  it  has  experienced. 
Thus,  Nt  >  0. 

9robability  of  entering  zeroth  (sixth)  cell,  having  had 
impacts  with  the  i*h  impact  occurring  in  cell  j. 

3robobility  of  (•)•  For  example,  P(j)  denotes  probo- 
tility  that  the  neutron  will  experience  exactly  j 

mpaets. 

’robability  of  baekseattering  with  the  last  impact 
n  the  j1*1  cell. 

)istance  (scaler)  measured  in  same  units  as  the 
'artesian  coordinates. 

Cnown  slab  thickness.  Thus,  space  occupied  by  the 
.lab  is  (-00  <  X  <  <»),  (-00  <  Y  <  ao),  (0  <  2  <  T)  . 

♦h  cell  entry  of  i*h  state  vector. 

-t  W  see  /3 . 
fb  state  vector. 

tn  arbitrarily  directed  one-dimensional  coordinate 
ystem  in  the  same  distance  units  as  the  Cartesian 
ystem,  and  with  origin  at  an  impact  point  of  the 
eutron. 

i  fixed  set  of  right-handed  Cartesian  axes  with  the 
litiol  slab  face  being  the  (X,  Y)  plane  and  the  -Z 
xis  being  totally  outside  the  slob. 

Inferential  albedo  per  steradian;  a  function  of 
‘-a,  2,,  T,  0O.  and  6,  where  the  integral  over  the 

ackward  (Z  <  0)  hemisphere  /nue  adft  =  (2  or, 

.  2rr  ,  -n/7  a 

quivalently,  Jq  Jq  Clsin  0  d0  dtp  =  U.  . 

'ifferentiol  scatter-through  (collided  transmission) 
er  steradian;  a  function  of  2„,  T,  P0,  and  0, 
here  the  integral  over  the  forward  (Z  >  T)  hemi- 
phere  /fhs  a’dft  =  (2  -  (?u  or,  equivalently, 

f  w  a*  sin  04044>  =  G*  -  • 

)  Jrr/1 

ngle  between  a  ray  and  the  positive  W  axis. 

ngle  between  projection  of  roy  onto  a  plane  perpen- 
ieular  to  the  W  axis  at  W  ~  0  ond  a  fixed  direction 
i  this  plane. 


0,  <t>  Polar  ond  azimuthal  angles  of  neutron  return  —  if  it 

does  return.  If  X*  Y!  Z‘  is  a  new  set  of  axes  (old 
set  translated  but  not  rotated),  having  its  origin  at 
the  neutron  return  point,  then  0  is  the  polar  angle 
that  the  returned  neutron  mokes  with  -Z  axis,  and 
is  the  azimuthal  angle  that  the  returned  neutron 
makes  with  -X’  axis;  where,  respectively,  the  direc¬ 
tions  -X1,  Y1,  X*,  ond  -Y*  have  azimuths  <t>  equal  to 
0,  rr/2,  w,  ond  3^/2.  When  0  -  0,  4>  is  undefined. 

0O  Polar  angle  of  neutron  incidence.  Neutron  is  directed 

from  outside  the  slob  without  interference  toward  the 
Cartesian  origin  (0,  0,  0).  This  direction  of  flight  is 
such  that  it  mokes  angle  0O  with  -Z  oxis,  angle 
(•77/2  -  0O)  with  -X  axis,  and  angle  rr/1  with  Y  axis; 

o  <  e0  <  77/2 . 

2*,  Z,,Xa  Known  constants;  respectively,  total,  scattering, 
and  absorption  macroscopic  cross  sections  of  slab 
material.  ■  Distance  traveled  by  the 

neutron  within  the  slab  between  two  successive 
impacts  is  a  random  variable  with  statistical  density 

function,  f(S)  =  A  priori  probability  of 

being  absorbed  at  the  next  impact  is  Sg  /^t  ■  The 
units  for  each  2  are  the  some  as  for  1/S  . 


<t>o  Azimuthal  angle  of  neutron  incidence;  by  definition, 

always  zero  except  when  0O  -  0,  and  then  <Po  i» 
undefined. 

n  Solid  angle  measured  in  steradians;  where,  e.g., 
ft  =  2t7  for  a  hemisphere. 
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YORK  09597 

COMMANDING  OFFICER,  ATTN  PUBLIC  WORKS  OFFICER.  U.  S.  NAVAL  STATION,  FPO  NEW 
YORK  09555 

COMMANDING  OFFICER,  U.  S.  NAVAL  STATION,  FPO  NEW  YORK  09571 

CHIEF  OF  ENGINEERS.  U.  S.  ARMY,  ATTN  ENGCW-E .  WASHINGTON  D.  C.  20315 

CHIEF  OF  ENGINEERS.  U.  S.  ARMY,  ATTN  ENGMC-E .  WASHINGTON  D.  C.  20315 


HEADQUARTERS.  U.  S.  AIR  FORCE.  DIRECTORATE  OF  CIVIL  ENGINEERING.  ATTN  AFOCE-ES. 
WASHINGTON  D.  C.  20330 

COMMANDING  OFFICER.  U.  S.  NAVAL  CONSTRUCTION  BATTALION  CENTER.  ATTN  MATERIEL 
DEPARTMENT,  CODE  140,  PORT  HUENEME,  CALIF.  93041 

DIRECTOR.  COAST  AND  GEODETIC  SURVEY,  U.  S.  DEPARTMENT  OF  COMMERCE,  6001 
EXECUTIVE  BOULEVARD.  ROCKVILLE.  MD.  20852 

DIRECTOR  OF  DEFENSE  RESEARCH  AND  ENGINEERING.  ROOM  3C-128,  THE  PENTAGON,  ATTN 
TECHNICAL  LIBRARY,  WASHINGTON  D.  C.  20301 

U.  S.  BUREAU  OF  RECLAMATION,  DEPARTMENT  OF  INTERIOR,  ATTN  MR.  T.  W.  MERMEL  , 
WASHINGTON  D.  C.  2C240 

FACILITIES  OFFICER,  ATTN  CODE  10B,  OFFICE  OF  NAVAL  RESEARCH,  WASHINGTON  D.  C. 
20360 

COMMANDER  NAVAL  BEACH  GROUP  TWO,  ATTN  PROJECT  OFF1CFR,  U.  S.  NAVAL  AMPHIBIOUS 
BASE.  LITTLE  CREEK,  NORFOLK,  V A.  23521 

U.  S.  ARMY,  ENGINEER  RESEARCH  AND  DEVELOPMENT  LABORATORIES.  ATTN  STINFO  BRANCH, 
FORT  BELVOIR,  VA.  22060 

AIR  FORCE  WEAPONS  LABORATORY,  KIPTLAND  AIR  FORCE  BASE.  ATTN  CODE  WLRC. 
ALBUQUERQUE.  N.  MCX.  87117 

LIBRARY,  DEPARTMENT  CF  METEOROLOGY  AND  OCEANOGRAPHY,  U.  S.  NAVA_  POSTGRADUATE 
SCHOOL.  .MONTEREY,  CALIF.  9394U 

LIBRARY  OF  CONGRESS.  WASHINGTON  D.  C.  2036C 

LIBRARY,  PUBLIC  DOCUMENTS  DEPARTMENT,  DUKE  UNIVERSITY,  DURHAM,  N.  C.  27706 

LIBRARY,  CIVIL  ENGINEERING  DEPARTMENT.  UNIVERSITY  OF  HAWAII,  HONOLULU,  HAWAII 
06822 

DIRECTOR,  INSTITUTE  OF  FISHERIES  RESEARCH,  UNIVERSITY  OF  NORTH  CAROLINA, 

MOREHEAD  CITY.  N.  C.  28557 

DIRECTOR.  AIR  UNIVERSITY  LIBRARY,  ATTN  AUL3T,  MAXWELL  AIR  FORCE  BASE,  ALA. 

36112 

OFFICER  IN  CHARGE,  CECOS,  ATTN  NDCE  COURSE.  PORT  HUENEME,  CALIF.  93041 

DEFENSE  INTELLIGENCE  AGENCY,  ATTN  DJAAP-1K2,  WASHINGTON  0.  C.  20301 

SPACE  SYSTEMS  DIVISION,  LOS  ANGELLS  AIR  FORCE  STATION,  ATTN  SSSD,  LOS  ANGELES. 
CALIF.  90C45 

COMMANDING  OFFICER.  'J.  S.  ARMY,  ENGINEER  RESEARCH  AND  DEVELOPMENT  LABORATORIES, 
ATTN  TECHNICAL  DOCUMENTS  CENTER,  FORT  BELVOIR.  VA.  22060 

DIRECTOR,  U.  S.  ARMY  ENGINEER  WATERWAYS  EXPERIMENT  STATION,  ATTN  MR,  G.  L. 
ARHUThNOT  JR,  P.  0.  bOX  631.  VICKSBURG,  miss.  39181 

U.  S.  ARMY  NUCLEAR  DEFENSE  LAoOPATORY,  ATTN  LIBRARIAN,  EDGEwOOD  ARSENAL,  MD. 

2  10  1 C 

U.  S.  ATOMIC  ENERGY  COMMISSION,  DIVISION  OF  TECHNICAL  INFORMATION,  P.  0.  BOX  62. 
OAK  RIDGE.  TENN.  37831 

CHIEF,  CIVIL  EFFECTS  BRANCH,  DIVISION  OF  BIOLOGY  AND  MEDICINE,  U.  S.  ATOMIC 


ENERGY  COMMISSION.  WASHINGTON.  0.  C.  20545 


DISASTER  RECOVERY  TRAINING  OFFICER,  U.  S.  NAVAL  CONSTRUCTION  BATTALION  CENTER. 
ATTN  CODE  TOC.  DAVISVILLE.  R.  I.  02854 

COMMANDING  OFFICEH.  BALLISTICS  RESEARCH  LABORATORIES,  ATTN  AMXBR-TB 
(MR.  WILLIAM  J.  TAYLOR).  ABERDEEN  PROVING  GROUND,  MD.  21005 

LCDR  CHARLES  W.  GULICK  JR.  CEC,  USN,  CNAbATRA  PUBLIC  WORKS.  NAS  PENSACOLA, 

FLA.  32508 

CAPT  W.  M.  MCLELLON,  CEC.  USN.  RET.,  468  1ST  STREET.  ALBANY,  N.  Y.  12210 

COR.  J.C.  LEDOUX,  NAVAL  FACILITIES  ENGINEERING  COMMAND,  WASHINGTON,  D.C.  20390 

LT  EDWARD  S.  PERRY.  U.S.  NAVAL  ORDNANCE  TEST  STATION,  3202  FOOTHILL  BLVD., 
PASADENA,  CALIF.  91107 

CAPT  L.N.  SAUNDERS.  CEC.  USN.  NAVAL  CIVIL  ENGINEERING  LABORATORY,  PORT  HUENEME. 
CALIF.  93041 

CDR  W.  A.  WALLS.  CEC.  USN,  U.  S.  NAVAL  MOBILE  CONSTRUCTION  BATTALION  FOUR 
FPO  SAN  FRANCISCO  96601 

LT  COL  C.  D.  DANIEL  (USA),  NUCLEAR  BRANCH  CHIEF,  DEFENSE  ATOMIC  SUPPORT  AGENCY. 
WASHINGTON  D.  C.  20310 

CDR  J.  F.  DOBSON,  CEC,  USN.  OFFICE  OF  CIVIL  DEFENSE.  DEPARTMENT  OF  THE  ARMY, 
WASHINGTON  D.  C.  20301 

CDR  J.  D.  ANDREWS.  CEC.  USN,  EXECUTIVE  OFFICER,  U.  S.  NAVAL  CIVIL  ENGINEERING 
LABORATORY,  PORT  HUENEME,  CALIF.  93041 

LCDR  GEORGE  M.  GANS.  CEC,  USN.  U.  S.  NAVAL  SCHOOL.  CEC  OFFICERS,  PORT  HUENEME, 
CALIF.  93041 

CDR  P.  A.  PHELPS.  CEC,  USN.  U.  S.  NAVAL  MCB  EIGHT,  C/0  FPO  SAN  FRANCISCO 
96601 

CDR  EDWARD  M,  SAUNDERS.  CEC.  USN,  NORTHWEST  DIVISION  NAVAL  FACILITIES 
ENGINEERING  COMMAND.  1633  W.  LAWTON  WAY,  SEATTLE.  WASHINGTON  98119 

LCDR  T.  YOSHIHARA,  CEC.  USN.  DOICC,  SOUTHEAST  ASIA,  APO  SAN  FRANCISCO.  96243 


LCDR  WALTLR  J.  EAGER  JR,  CEC,  USN.  POSITIVE  ION  ACCELERATOR  BRANCH,  AFRR  J -NNMC  , 
BETHESOA,  MD.  20C14 

LT  L.  K.  DONOVAN,  CEC,  USN,  DISASTER  CONTROL  OFFICER,  U.  S.  NAVAL  COMMUNICATION 
STATION,  HONOLULU.  FPC  SAN  FRANCISCO,  CALIF.  96613 

CDR  W.  J.  FRANCY,  CEC.  USN.  ENGINEERING  DEPARTMENT,  U.  S.  NAVAL  ACADEMY, 
ANNAPOLIS.  «D.  21402 

LCDR  1.  0 •  CROWLEY  JP,  BLAST  AND  SHOCK  DIVISION.  DEFENSE  ATOMIC  SUPPORT 
AGENCY,  WASHINGTON  0.  C.  2031C 

LCDR  w.  H.  BANNISTER,  CEC,  USN,  FIELD  COMMAND,  DEFENSE  SUPPORT  AGENCY.  SANDIA 
BASE.  ALbUOUEROUE.  N.  MEX.  87115 

MAJOR  RObERT  E.  CRAWFORD.  AIR  FORCE  WEAPONS  LABORATORY  (WLDC).  K I PTL  AND  AIR 
FORCE  BASE.  N.  MEX.  87117 


CDR  C.  R.  WHIPPLE.  CEC.  USN ,  U.  S.  NAVY  PUBLIC  WORKS  CENTER.  BOX  6.  FPO  SAN 
FRANCISCO.  CALIF.  96630 

COR  ALLEN  F.  DILL.  CEC.  USNR.  FIELD  COMMAND,  DASA  SANDIA  BASE 
ALBUOUERQUE,  NEW  MEXICO  87115 

MR.  NEAL  FIT2SIM0NS  IOISECTORI.  ENGINEERING  DEVELOPMENT  DIVISION,  OFFICE  OF 
CIVIL  DEFENSE,  DEPARTMENT  OF  DEFENSE,  WASHINGTON  D.  C.  20310 

MR.  BENJAMIN  C.  TAYLOR.  SPECIAL  ASSISTANT  TECHNICAL  LIAISON,  OFFICE  OF  CIVIL 
DEFENSE,  DEPARTMENT  OF  ARMY,  WASHINGTON  D.  C.  20310 

MR.  CHARLES  M.  EISENHAUER,  RADIATION  PHYSICS  LABORATORY.  NATIONAL  BUREAU  OF 
STANDARDS.  WASHINGTON  D.  C.  20234 

DR.  JAMES  0.  BUCHANAN  (DIRECTOR).  SHELTER  RESEARCH,  OFFICE  OF  CIVIL  DEFENSE. 
DEPARTMENT  OF  DEFENSE.  WASHINGTON  D.  C.  20310 

MR.  JOHN  A.  AUXIER.  OAK  RIDGE  NATIONAL  LABORATORY 

MR.  GEORGE  N.  SISSON,  OFFICE  OF  CIVIL  DEFENSE.  DE 
WASHINGTON  D.  C.  20310 

MR.  F.  J.  TAMANINI,  A  AND  E  StRVICES  DIVISION,  OFFICE  OF  CIVIL  DEFENSE. 
DEPARTMENT  OF  DEFENSE.  WASHINGTON  D.  C.  20310 

MR.  CHARLES  E.  CLIFFORD,  OAK  RIDGE  NATIONAL  LABORATORY,  P.  0.  BOX  X 
OAK  RIDGE.  TENNESSEE  37831 

MR.  WILLIAM  R.  PERRET,  SANDIA  LABORATORY,  DIVISION  5232,  SANDIA  CORPORATION, 
ALBUQUERUUE,  N.  MEX.  87116 

MR.  RICHARD  PARK,  ADVISORY  COMMITTEE  ON  CIVIL  DEFENSE.  NATIONAL  ACADEMY  OF 
SCIENCES.  2101  CONSTITUTION  AVENUE  N.  W.,  WASHINGTON  D.  C.  20418 

DR.  E.  E.  MASSEY,  DEFENSE  RESEARCH  BOARD,  OFFICE  OF  DIRECTOR  OF  ATOMIC  RESEARCH, 
DEPARTMENT  OF  NATIONAL  DEFENSE.  OTTAWA,  CANADA 

DR.  NICHOLAS  PERRONt,  OFFICE  OF  NAVAL  RESEARCH,  ATTN  CODE  439,  DEPARTMENT  OF 
THE  NAVY.  WASHINGTON  D.  C.  20390 

MR.  J.  G.  JACKSON  JR  (CHIEF),  IMPULSE  LOADS  SECTION.  SOIL  DYNAMICS  BRANCH, 

U.  S.  ARMY  ENGINEER  WATERWAYS  EXPERIMENT  STATION,  VICKSBURG,  MISS.  39180 

i 

LIBRARY,  PROTECTIVE  STRUCTURES  DEVELOPMENT  CENTER,  BUILDING  2591.  FORT  BELVOIR. 
VA.  22060 

LIBRARIAN,  NATIONAL  RESOURCE  EVALUATION  CENTER,  OFFICE  OF  EMERGENCY  PLANNING. 
EXECUTIVE  OFFICE  BUILDING,  WASHINGTON  D.  C.  20504 

MR.  ERWIN  E.  SHALOWITZ,  PROTECTIVE  CONSTRUCTION.  GSA  BUILDING,  19TH  AND  F 
STREETS  N.  W.,  WASHINGTON  0.  C.  204C5 

DR.  WILLIAN  E.  KREGER  (HEAD).  NUCLEONICS  DIVISION.  U.  S.  NAVAL  RADIOLOGICAL 
DEFENSE  LABORATORY,  SAN  FRANCISCO,  CALIF.  94135 

DR.  WILLIAM  L.  WHITE,  DIRECTOR,  CIVIL  DEFENSE  TECHNICAL  OFFICE.  STANFORD 
RESEARCH  INSTI1UTE,  MENLO  PARK,  CALIF.  94C25 

CHIEF  ENGINEERING  DIVISION.  U.  S.  ARMY  ENGINEER  DISTRICT  (OMAHA),  CORPS  OF 
ENGINEERS,  6012  U.  S.  POST  OFFICE  AND  COURT  HOUSE,  OMAHA,  NEB.  68102 

MR.  FRANK  KOPTONAK,  U.  S.  NAVY  UNDERWATER  SOJND  LABORATORY,  NEW  LONDON,  CONN. 


.  OAK  RIDGE.  TENN.  37830 

PARTMENT  OF  DEFENSE, 

h  i 


06320 


COMMANDER.  U.  S.  NAVAL  WEAPONS  LABORATORY,  ATTN  CODE  TEV,  DAHLGREN,  VA.  22448 

MR.  DAVID  A.  PALMER.  CHESAPEAKE  DIVISION.  NAVAL  FACILITIES  ENGINEERING  COMMAND. 
WASHINGTON  NAVY  YARD.  WASHINGTON.  D.  C.  20390 

LIBRARY,  DEPARTMENT  OF  METEOROLOGY  AND  OCEANOGRAPHY.  U.  S.  NAVAL  POSTGRADUATE 
SCHOOL,  MONTEREY,  CALIF.  93940 

DR.  C.  V.  CHELAPATI,  ASSISTANT  PROFESSOR  OF  ENGINEERING,  CALIFORNIA  STATE 
COLLEGE  AT  LOS  ANGELES,  LOS  ANGELES,  CALIF.  9C032 

MR.  G.  H.  ALBRIGHT,  DEPARTMENT  OF  ARCHITECTURAL  ENGINEERING.  101  ENGINEERING 
•A'  BUILDING,  THE  PENNSYLVANIA  STATE  UNIVERSITY,  UNIVERSITY  PARK,  PA.  16802 

PROFESSOR  J.  NEILS  THOMPSON,  EDUCATION  AND  RESEARCH,  CIVIL  ENGINEERING 
DEPARTMENT,  UNIVERSITY  OF  TEXAS.  AUSTIN,  TEX.  78712 

MR.  FRED  M.  SAUER,  PHYSICS  DEPARTMENT,  STANFORD  RESEARCH  INSTITUTE.  MENLO  PARK, 
CALIF.  940 2 S 

DR.  THEODORE  SCHIFFMAN.  ILLINOIS  INSTITUTE  OF  TECHNOLOGY.  RESEARCH  INSTITUTE. 

10  W.  35Th  STREET.  CHICAGO,  ILL.  60616 

DR.  ROBERT  V.  WHITMAN,  ROOM  1-346.  MASSACHUSETTS  INSTITUTE  OF  TECHNOLOGY, 
CAMBRIDGE.  MASS.  02139 

DR.  LEWIS  V.  SPENCER,  PHYSICS  DEPARTMENT,  OTTAWA  UNIVERSITY,  OTTAWA,  KANSAS 
66067 

MR.  WERNER  WEBER,  DIRECTOR  OF  RESEARCH  AND  DEVELOPMENT,  PUBLIC  SECURITY 
BUILDING.  STATE  OFFICE  BUILDING  CAMPUS.  ALbANY,  N.  Y.  12226 

MR.  RICHARD  D.  CAVANAUGH,  3APRY  CONTROLS  INC..  70C  PLEASANT  STREET,  WATERTOWN 
MASS.  02172 

MR.  KENNETH  KAPLAN,  URS  CORPORATION,  1811  TROUSDALE  DRIVE,  BURLINGAME,  CALIF. 
94011 

MR.  WALTER  D.  GUNTHER,  GOVERNMENT  CONSULTING,  THE  MITRE  CORPORATION,  P.  0.  BOX 
208,  LEXINGTON,  MASS.  02173 

MR.  LYNDON  WELCH,  EBERLE  M.  SMITH  ASSOCIATES  INC.,  153  EAST  ELIZABETH  STREET, 
DETROIT,  MICH.  482C1 

DR.  MERIT  P.  WHITE  (HEAD).  DEPARTMENT  OF  CIVIL  ENGINEERING,  UNIVERSITY  OF 
MASSACHUSETTS,  AMHERST.  MASS.  01003 

DR.  ROBERT  J.  HANSEN,  DEPARTMENT  OF  CIVIL  AND  SANITARY  ENGINEERING, 
MASSACHUSETTS  INSTITUTE  OF  TECHNOLOGY,  CAMfcRIDC-E,  MASS.  02138 

MR.  LUKE  J.  VORTMA.N,  DIVISION  5232  ,  SAND  I A  CORPORATION,  SANDIA  bASE, 
ALBUQUERQUE,  N.  MEX.  87116 

DR.  ERIC  T.  CLARKE.  VICE  PRESIDENT  FOR  RESEARCH  AND  DEVELOPMENT,  TECHNICAL 
OPERATIONS  INC..  BURLINGTON,  MASS.  02103 

DR.  A.  B.  CHILTON  (ASSOCIATE  PROFESSOR  I  ,  CIVIL  AND  NUCLEAR  ENGINEERING,  207 
NUCLEAR  ENGINEERING  LABORATORY,  UNIVERSITY  OF  ILLINOIS,  URBANA ,  ILL.  618C3 

MRS.  SHEA  VALLEY,  ERTE .  A.  F.  CAMBRIDGE  RESEARCH  LABORATORIES,  L.  G.  HANSCOM 
FIELD.  BEDFORD.  MASS.  01731 


DR.  JOHN  r.  HANLEY.  ASSOCIATE  PROFESSOR  OF  CIVIL  ENGINEERING*  UNIVERSITY  OF 
MINNESOTA.  MINNEAPOLIS.  MINN.  55414 

PROFESSOR  JOSEPH  SILVERMAN,  DEPARTMENT  OF  CHEMICAL  ENGINEERING,  UNIVERSITY  OF 
MARYLAND,  COLLEGE  ftARK,  MD.  20742 

OR.  F.  T.  MAVIS.  DEAN  AND  PROFESSOR,  DEPARTMENT  OF  MECHANICAL  ENGINEERING, 
UNIVERSITY  OF  MARYLAND,  COLLEGE  PARK,  MO.  20742 

PROFESSOR  MATEO  L.  P.  GO,  PROFESSOR  OF  CIVIL  ENGINEERING,  DEPARTMENT  OF  CIVIL 
ENGINEERING,  UNIVERSITY  OF  HAWAII,  HONOLULU,  HAWAII  96822 

DR.  JAMES  P.  ROMUALDI,  PROFESSOR,  RESEARCH  IN  STRUCTURAL  DYNAMICS.  CARNEGIE 
INSTITUTE  OF  TECHNOLOGY.  PITTSBURG,  PA.  15213 

DR.  ALEKSANDAR  SEDMAK  VESIC,  PROFESSOR  OF  CIVIL  ENGINEERING,  DUKE  UNIVERSITY, 
DURHAM,  N.  C.  27706 

MR.  C.  C.  MOW,  THE  RAND  CORPORATION,  1700  MAIN  STREET.  SANTA  MONICA,  CALIF. 
90406 

DR.  HAROLD  L.  BROOfc.  THE  RAND  LORPORAT ION, 1700  MAIN  STREET.  SANTA  MONICA, 
CALIF.  90406 

PROFESSOR  RAYMOND  R.  FOX.  DIRECTOR,  NUCLEAR  DEFENSE  ENGINEERING  COURSES.  THE 
GEORGE  WASHINGTON  UNIVERSITY,  WASHINGTON  D.  C.  20006 

DR.  S.  H.  LEVINE  (HEAD),  NUCLEAR  SCIENCES  LABORATORY,  NORTHROP  SPACE 
LABORATORIES,  3401  WEST  BROADWAY,  HAWTHORNE,  CALIF.  90250 

PROFESSOR  ELMIRA  SMYRL ,  SCHOOL  OF  ARCHITECTURE,  MONTANA  STATE  COLLEGE. 
BOZEMAN,  MONT.  59715 

PROFESSOR  CARL  H.  KOONPV.  CIVIL  ENGINEERING  DEPARTMENT,  WORCESTER  POLYTECHNIC 
INSTITUTE.  WORCESTER.  MASS.  01609 

DR.  ARMAS  LAUPA,  THE  RAND  CORP,  1700  MAIN  STREET,  SANTA  MONICA,  CALIF.  90406 

DR.  C.  P.  S I  ESS ,  DEPARTMENT  OF  CIVIL  ENGINEERING.  UNIVERSITY  OF  ILLINOIS, 
URBANA,  ILL.  61803 

PROFESSOR  ADRIAN  F.  RICHARDS,  GEOLOGY  AND  CIVIL  ENGINEERING,  DEPARTMENT  CF 
GEOLOGY,  UNIVERSITY  OF  ILLINOIS,  URBANA,  ILL.  61803 

PROFESSOR  W.  J.  HALL,  DEPT.  OF  CIVIL  ENGINEERING.  Ill  TALBOT  LABORATORY, 
UNIVERSITY  OF  ILLINOIS.  URBANA,  ILL.  61803 

AEROSPACE  CORPORATION,  SAN  BERNARD  I  NO  OPERATIONS.  P.O.  BOX  1308,  SAN 
BERNARDINO.  CALIF.  92402 

MR.  R.  M.  WEBB.  HOUSING  AND  HOME  FINANCE  AGENCY,  COMMUNITY  FACILITIES 
ADMINISTRATION.  WASHINGTON  D.  C.  20410 

DIRECTOR.  EASTERN  DIVISION,  BUREAU  OF  YARDS  AND  DOCKS,  90  CHURCH  STREET, 

NEW  YORK,  NEW  YORK  100C7 

CHIEF,  SPECIFICATIONS  BRANCH,  DESIGN  DIVISION,  PUBLIC  BUILDINGS  SERVICE  PCCX, 
GENERAL  SERVICES  ADMINISTRATION,  WASHINGTON,  D.  C.  20405 

COMMANDANT,  U.S.  ARMY  LOGISTICS  MANAGEMENT  CENTER.  FORT  LEE,  VIRGINIA  23801 
ATTN  DEFENSE  LOGISTICS  STUDIES  INFORMATION  EXCHANGE 

COMMANDING  OFFICER,  DEPARTMENT  OF  THE  ARMY,  MISSOURI  RIVER  DIVISION, 


CC. Cl  •  .OU.Lt.RS,  P.O.  o'JX  103.  DOWNTOWN  STATION.  OMAHA.  NEBRASKA  68101 

Ar.RL  -  RESEARCH  LIBRARY,  ATTN  CRMXLR/STOP  29,  L.  G.  HANSCOM  FIELD,  BEDFORD, 
MASS.  01731 

ARCTIC  At  ROMEO  I  CAL  LABORATORY,  U.  S.  AIR  FORCE.  ATTN  LIBRARY,  APO  SEATTLE. 
WASH.  98731 

COMMANDING  OFFICER.  U.  S.  A PMY  BIOLOGICAL  LABORATORIES.  ATTN  DOCUMENTS  - 
TECHNICAL  LIBRARY,  FORT  DFTcjCK,  FREDERICK,  MD.  21701 

■MR.  P.  C.  SMITH,  PRINCIPAL  RESEARCH  ENGINEER  (SOILS).  MATERIALS  DIVISION. 

BUREAU  OF  PUBLIC  ROADS.  WASHINGTON  D.  C.  20235 

MR.  Ju-.N  H.  HUoNLLL.  RADIATION  PHYSICS  DIVISION,  NATIONAL  BUREAU  OF  STANDARDS. 
WASHINGTON  0.  C.  20234 

MR.  S.  K.  PENNY  (DIRECTOR),  RADIATION  SHIELDING  INFORMATION  CENTER,  OAK  RIDGE 
NATIONAL  LABORATORY,  P.  0.  BOX  X.  OAK  RIDGE,  TENN.  37831 

officer  in  charge',  u.  s.  naval  biological  laboratory,  naval  supply  center, 

OAKLAND,  CALIF.  94625 

COMMANDING  OFFICER,  NRCTC  UNIT,  RENSSELAER  POLYTECHNIC  INSTITUTE.  TROY.  N,  Y, 
12181 

CHIEF.  tsL  (CAO  OF  MEDICINE  AY'  SURGERY.  ATTN  RESEARCH  DIVISION,  NAVY  DEPARTMENT 
WASHINGTON  D.  C.  20390 

OFFICER  IN  CHARGE,  U.  S.  NAVAL  SUPPLY  RESEARCH  AND  DEVELOPMENT  FACILITY.  ATTN 
LIBRARY,  NAVAL  SUPPLY  CENTS  R  ,  SAVONNE.  N,  J.  07002 

CHIEF  CHEMIST,  ATTN  CODE  32.N5,  NORFOLK  NAVAL  SHIPYARD,  PORTSMOUTH,  VA.  23709 

MR.  ROSS  E.  MORRIS,  RUBBER  LABORATORY,  MARE  ISLAND  DIVISION,  SAN  FRANCISCO  BAY 
NAVAL  SHIPYARD,  SAN  FRANCIS'S.  CAL  IF,  94135 

MR.  W.  A.  2ISMAN,  CHEMISTRY  [VISION.  ATTN  CODE  6100,  U.  S,  NAVAL  RESEARCH 
LABORATORY ,  WASHINGTON  0.  C.  20390 

LlePAFY,  J.  S.  ARMY,  COLD  REGIONS  PESEARCH  AND  ENGINEERING  LABORATORY,  P.  0. 

BOX  282.  hANCV'R,  N.  H.  03755 

MR.  K.  .  WlL  'UR.  PROFESSOR  CP  ZOOLOGY,  OUKE  UNIVERSITY.  DURHAM,  N.  C.  27706 

DR.  CHARLES  :.  LANE,  PROFESSOR.  INSTITUTE  OF  MARINE  SCIENCE.  UNIVERSITY  OF 
MIAMI,  CORAL  CABLES.  FLA.  33146 

LIBRARY,  DOCUMENTS  DIVISION,  PRINCETON  UNIVERSITY,  PRINCETON,  N.  J,  08540 
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